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Fields from a relativistic magnetic explosion 



00 

o 
o 

(N 
D 



O 



> 

m 
o 

G\ 
O 
oo 
O 



K.N. Gourgouliatos* & D. Lynden-Bellf 

Institute of Astronomy, University of Cambridge, Madingley Road CBS OHA 



Accepted 2008 August 25. Received 2008 August 25; in original form 2008 July 16 



ABSTRACT 

Following Prendergast we study the relativistically expanding electromagnetic fields 
generated by an axisymmetric explosion of magnetic energy in a small volume. The 
magnetic field expands uniformly either within a cone or in all directions and it is 
therefore accompanied by an electric field. In the highly conducting plasma the charges 
move to annul the electric field in the frame of the moving plasma. The solutions pre- 
sented are analytical and semi-analytical. We find that the time-scale for the winding 
up of the initial magnetic field is crucial, as short time-scales lead to strong radiant 
fields. Assuming a magnetic field of lO^^Gauss emerging from a magnetosphere of 
lO^cm we end with a jet when confined by a pressure environment that falls more 
slowly than r~^. The jet carries energy of lO^^erg, which is mostly due to differential 
rotation at the base. 
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1 INTRODUCTION 

Observations of a wide variety of astronomical objects 
suggests the existence of magnetic fields in relativistic 
environments. There have been many studies of magnetic 
fields emanating from differentially rotating systems. Some 
of them confine themselves to the non-relativistic regime 
in which the displacement currents ca n be ne g lected . 
Even then few of th e m are a nalytic e.g. iLovelacd 
Lvnden-Bell fc Boilvl (|l994l ). ISautv fc Tsinganod 



197e), 



19941 ) ■ 



Lvnden-Belll (|2006|). while most of them are computa- 



tional e.g . 'Le Blanc & Wilson' ('l970l). iBell fc Luce^ l|l995h . 
lOuved. Pudritz fc Stone (|1997, ). In the force- free case it 
was shown that the time evolution of the magnetic field 
arises solely from the time dependence of the boundary 
conditions so that the exact dynamical evolution can be 
calculat ed from the time dependent sequence of static 
models (|Lvnden-Bel]| [2009) . Those in turn can be derived 
from the energy principle. However, that simplification 
depends on both the force-free condition and the neglect 
of the displacement currents which is only valid when the 
velocities are much less than c. Relativistic problems are 
harder as that approximation is i nvalid, so most studies 
pu rely comput a tional e.g. Li. Chiueh fc Begelman 



2002) 



1992|) . iKomissarov' ("2002'), 'Gammie. McKinnev fc Toth 



e Villicrs, Hawlcy, Krolik & Hiroso (2003" 

l|2008l ) or f 



iTchekhovskov. McKinnev fc NaravanI 
analytical e.g. Contopoulo j ( 19951 ) 



and 



* E-mail: kng22@ast.cam.ac.uk 
t E-mail: dlb@ast.cam.ac.uk 



Despite this, iPrendergastI (|2005l ) was able to find an 
exact solution to the relativistic MHD problem of a point 
magnetic explosion. He derived a time dependent relativis- 
tic analogue of the Grad-Shafranov equation that governs 
axially symmetric force-free MHD by assuming that the ra- 
dial coordinate and the time since the explosion only appear 
in the dimensionless combination v = r/ct. Strictly speak- 
ing such equations are only valid when the length scale and 
cx the time scale of the region of the explosion are much 
smaller than r and ct. The resulting Prendergast equation is 
non-linear but becomes linear in a special case. It was this 
special case that Prendergast studied in detail. However, he 
found that there are spherical nodes where the radial mag- 
netic field is zero and these nodes occur before the highly 
relativistic regime r « ct is reached. This has the unfortu- 
nate consequence that magnetic field lines emanating from 
small r turn back before they reach the extremely relativistic 
region where the displacement currents are very important. 
The field lines in that region are an appendage unattached 
to their origin. Very similar effects are well known when one 
takes the analogous case, j = aB with a constant, in non- 
relativistic MHD in spherical coordinates. When ar becomes 
large, the system gives way to oscillating solutions with a se- 
ries of nodes. This difficulty occurs because has dimen- 
sion L and the field has to vary on this fixed length scale 
even at large r. At large distances there is too much current 
for unit field and the smaller scale is then reflected in the 
scale of oscillation. 

It is known that non-linear ansatzes for the non- 
relativistic Grad-Shafranov equation can avoid this, which 
we now recognise as a bad consequence of a mathematically 
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simple linear approximation which is not generally justified 
in the physics of the problem. We shall therefore study the 
non-linear Prendergast equation in all its glory! 

We generalise the idea of force-free fields by considering 
a configuration where for each point there exists a Lorentz 
frame in which the electric field is zero and the current is 
along the magnetic field. The force density f in the frame 
fixed at the origin is then: 



f = p(E + V X B) = 0, 



(1) 



where cv(r, t) is the velocity of the moving frame. In 
this pape r, we study such fie lds by solving the equation pro- 
posed by IPrendergastI (|2005l ). However, before stepping to a 
solution, we extract as much information as possible about 
the fields that is independent of the detailed form of the so- 
lution. Then we solve the equations semi-analytically as it 
impossible to achieve general analytical solutions. The solu- 
tions si mplify in the non-rela t ivistic limit and converge to 
those of iLvnden-Bell fc Boilvl (|l994h . There are some other 
cases that are interesting and exactly soluble, namely the 
current-free magnetic dipole and the linear force-free field 
of Prendergast. The introduction of currents in the system 
allows the existence of a toroidal component of the field at 
the cost of making analytical solutions much harder, how- 
ever it is still possible to design analytical solutions for this 
structure. 



2 PROBLEM SETUP AND SOLUTION 
STRATEGY 

2.1 Self-Similar form of Maxwell's Equations 

In accordance to Prendergast's formalism we consider an 
electromagnetic field configuration for which there is at least 
one frame of reference at each point of space where the 
electric field vanishes and the magnetic field is parallel to 
the electric current. We denote these frames and the phys- 
ical quantities appearing there by a prime. We take each 
of these frames to move uniformly with velocity v = r/ct. 
Although each frame moves uniformly away from the origin 
and at each moment those frames further from the origin 
move faster, we are not using the expanding coordinates of 
cosmology. The frame of reference of the observer at the ori- 
gin is the unprimed frame. The primed frames of reference 
move radially with respect to the unprimed frame under the 
scaled velocity v 



V — — r. 

ct 



(2) 



Then we may write for the primed frame of reference: 



aB', 



(3) 



E' = 0. (4) 
The magnetic field in the unprimed frame of reference 



Br ~ B'r, 



■yB'g, B4, = 75^, 



(5) 



where 7 = (1 — v"^) The electric field in the un- 
primed frame of reference is 

Er =0, Eg^ vB^, E^ = --vBg. (6) 
The electric current in the unprimed frame of reference 



jr = ^j'r + VCp, jg = j'g , = j'^. (7) 

By virtue ofV-B = 0, VxE = -dB/d{ct), the self- 
similarity condition and (1) we re-express the fields in terms 
of the fiux function P and T. P{r, 9, t) is the magnetic fiux 
through the cap of a sphere of radius r subtending a semi- 
angle 6 at the origin, and T is related to the toroidal com- 
ponent B^ and to the specific torque carried by field lines 
between P and P + dP. Remembering that v — r/ct and 
H = cos 6 we have 



B 



E 



dP_^ _ 
djj. 

vT 



,2 d 



T 



-.e 



dP - 



2 dv 



(8) 



(9) 



Now we substitute the currents (7) and the electric and 
magnetic fields into V x B = dE/d{ct) + 47r/cj and finally 
take into account charge conservation to get Prendergast's 
equation for the flux function 



dp ( v/3 



dp /■ v(i \ 



(10) 



In Prendergast's notation Ti.{P) is ^^(5{P) and our P 
and T are 2n times his. 

The toroidal part T is related to the source function /3 

by 



r = 



1 — 



(11) 



Our task is to determine the function P by solving the 
partial differential equation (10). This is not a straightfor- 
ward task; even under a numerical treatment of the problem, 
the solution depends on the form of the source function (5 
chosen. Our strategy for solving equation (10) is 

1. to study the non-relativistic limit so as to see how it con- 
nects to earlier work. 

2. to study the elementary case when there are no currents 
in the body of the plasma. The resulting self-similar solu- 
tions of Maxwell's equations become singular at r=ct if the 
fields at the origin are turned on suddenly. We demonstrate 
how these singularities disappear when the fields are turned 
on over a small finite time. 

3. find what form of source function I3{P) allows solutions 
of (10) that are separable in and v. The special form of 
/3(P) required can not be made to vanish at P = 0. How- 
ever it gives an exact solution in regions where P is above a 
certain threshold, Po, where /3(Po) = = /3'(Po). When P 
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is below that threshold the field joins smoothly to one with 
no currents. 

4. The current-free solutions in this outer cocoon are not of 
product form but are found by analysis into suitable har- 
monics whose coefficients are computed to high accuracy. 

5. Thus we end with a global solution that expands with the 
velocity of light within a narrow prescribed cone. It has high 
fields but no charges near r — ct. Below that we enter the 
current free cocoon, soon followed by the basic relativistic 
solution of product form which is accompanied by charges 
and currents. 



2.2 Properties of the field 

In this section we focus on the properties of the field. It may 
seem more natural to achieve a solution for equation (10) 
and then investigate the details of the field. However it is 
possible to obtain some physical understanding by studying 
general properties of the field that do not depend on the 
details of the solution. 

We seek frames of reference where the electric field van- 
ishes. Such frames are particularly interesting as they deter- 
mine the velocity of the field and allo w us to trace th e time 
evolution of the magnetic field lines jSturrockl 1 1994) . The 
setup of the problem gives an initial set of frames obeying 
this property, the uniformly expanding ones. Indeed the mo- 
tion of the field lines may be viewed as a simple expansion. 
A charge q that is placed at distance ro from the origin at 
time to with velocity vo = ro/ctor, does not feel any force, as 
there is no electric field and the magnetic field in this frame 
has no effect on a particle at rest in it. The particle will re- 
main on the same magnetic field line as those are determined 
by P which in turn is a function of v = r /ct. However, as the 
fields have components in the all three directions the velocity 
field of uniform expansion is not perpendicular to the mag- 
netic field. Conventionally the velocity of field lines is chosen 
to be perpendicular to them. We examine other frames of 
reference where there is no electric field. In order to find 
such frames we use the following quadratic equation for the 
velocity field cvf (|Landau fc Lifshitd [19751 ) ■ that describes 
frames of reference where the magnetic and the electric fields 
are parallel to each other 



VF 



l-hVF^ 



E X B 



(12) 



This velocity field is clearly perpendicular to B. Given 
that the electric and the magnetic field by construction sat- 
isfy E ■ B = 0, which is a relativistic field invariant, we 
conclude that in a frame where E' || B' either the electric 
or the magnetic field vanishes. Thus equation (12) may be 
used to define the velocity of the magnetic field. However 
the general velocity vf of a frame in which E' and B' are 
parallel allows boosts in the direction of the magnetic field 
and only determines the motion in the direction transverse 
to the magnetic field. The solution of the above quadratic 
equation gives the velocity we are interested in, however it 
is feasible to achieve a simpler formula for the motion of 
the field lines. Equation (1) yields that E = — v x B. Any 
velocity field of the form 



E X B 



Vf = 



IB! 



AB, 



(13) 



for arbitrary A, gives a set of frames of reference where 
the electric field vanishes . The Lorentz t ransformation of the 
electric field is given by (IJacksonlllQTSh 



E' = 7(E-f VF X B) 



7+1 



vf(vf ■ E) 



(14) 



Substituting (13) into (14) and working the cross prod- 
ucts we obtain that E' = 0. Therefore equation (13) gives 
the gener al velocities that can be att ributed to the field lines; 
we follow iLandau fc Lifshitd l|l975h in using the expression 
for A = which gives a velocity perpendicular to the mag- 
netic field. This is a solution to (12) provided |E| < |B|. 
By virtue of equations (8) and (9) we find that the usual 
velocity field is 



VF,r 



VF,9 = 



VF,4, — 



(dP_y 



I — ^"^ \ dv J 



+ v- 



dP dP 



47r2r4|BP 



1 



|BP ^1 - ^2 dfi dv 



vT dP 



47^2^41312 



(15) 



(16) 



(17) 



Relations (15)-(17) give the details of the motion of the 
field lines. At the first instance we can make some reasonable 
assumptions for the flux function P. A magnetic field con- 
sists of fleld lines rising from and returning to an imaginary 
spherical cap of radius ro and confined between angles and 
corresponding to fi = 1 and fi = fj.o.ln the case of a dipole 
field instead of a spherical cap we have the whole spherical 
surface, and no = —1. The flux function P is zero at the 
boundaries and has a single maximum at n — fii where the 
fleld lines turn back. The toroidal component T has to obey 
r — > approaching the edges and normally the decrease has 
to be such so there is no singularity. The above assumption 
allows the determination of the signs of the derivatives ap- 
pearing at relations (15)-(17). Therefore ^ 0, as there 
cannot be any flux generated as the field moves upwards; 
|£ > for ^0 M < Ati; If > for ^ii < sC 1 and |£ = 
at fi = fii. 

We conclude that the field lines viewed this way perform 
a complex motion that consists of a non-uniform radial ex- 
pansion; a meridional motion where the parts of the field 
lines that lie in /i € (/xi,l) move in the direction of —0, 
those in G (/io, Pi) in the direction of 0; and an azimuthal 
motion where the parts of the field lines lying in /i G (/lo, ^i) 
move in the direction of (f> and the rest in the direction of 
—(f). The magnitude of the velocity depends on the details 
of the solution and cannot be determined by the previous 
arguments. 

2.3 The non-relativistic limit 

The parameter v, that is the expansion velocity scaled to 
the speed of light gives an estimate of the importance of the 
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relativistic terms for the problem. In the hmit of u << 1 
equation (10) becomes, 



+ (1 - 



(18) 



Equation (18) is similar to the Grad-Shafranov equa- 
tion for axisymm etric magnetic fields (|Grad fc Rubinll 19581 : 
IShafranovlll966l '). There are man y solut i ons o f this equa- 
tion in astrophysical con text, IjPriestl Il984l : lAlvl 1 1994 
iLvnden-Bell fc Boilvl 1 19941 ). Amongst these solutions the 
most appropriate to our case is that of lLvnden-Bell fc Boilvl 
j 19941 ) where the field is a self-similar quadruple. In our first 
paradigm we are going to apply this solution to a self-similar 
dipole. 

Assuming a self-similar solution we set P = 
FmaxV~'' f{ii), where Fmax is a flux normalisation, formally 
the maximum fiux at u = 1, so that the angular part of the 
flux f{fi) < 1, therefore equation (18) becomes 



l{l + 1)V~' Fmax fill) + (1 - tl'')F^axV~' f{^J)" = (19) 

In order to achieve self-similar solutions we need 
fidp/dP oc ti~'~^. With P a product of a function of v and 
a function of ^, I3{P) can only be proportional to a power 
oivHp itself is a power of P, since P is proportional to v^^ 
the required power is given by 



P = c,P'+^'\ 

a more detailed derivation 
iLvnden-Bell fc Boilvl (|l994h . Setting co = 
(19) reduces to 



(20) 



ciFmix, equation 



/(/ + 1)/(m) + (1 - /i')/(M)" = + \)f (21) 

To solve equation (21) in a sphere we apply the bound- 
ary conditions /(I) = /(—I) = 0, and determine co by 
the condition that / has a single maximum f{fj.max) = 1- 
These conditions are sufficient to solve the equation for 
any I. For / = 1 the solution is /(/i) — 1 — fi^ and 
Co = 0, for I small the solution of (21) takes the form 
/ = 1 — 1/(1 + l)lncosh((l -I- Z~^)cop) and in general for 
intermediate values it is solved numerically. Therefore the 
fields take the following form in this non-relativistic limit 



E = 



27rr^i)' 

Frr 



/'f + 



If 



--0 + 



Co/ 



1 + 1/! 



(22) 



(23) 



The motion of the field lines can be found by equations 
(15)-(17) by applying the fields given at (22) and (23). 



Vf 



(1 - ^2)/'2 4. /2 J2 c2/2/'+2 ' 

+ (1 - fi^f^^lff'e + (1 - fiY^\f'f+'^'4,]. (24) 

This case is different to the Grad-Shafranov equation, as 



in the latter the field only depends on r whereas here the field 
depends on r/ct. However, a snapshot of the configuration 
for any given t satisfies the Grad-Shafranov equation. 



3 RELATIVISTIC SOLUTIONS 
3.1 Current-Free Solutions 

When there are no source terms on the right, equation (10) 
gives simple relativistic solutions with no toroidal magnetic 
field component. 

In this section we study this case; we show that one 
solution reduces to a linearly increasing magnetic dipole, 
where the field is zero outside a sphere expanding at the 
speed of light. 

With /3 = equation (10) reduces to the simple form 



By using Ogilvie's transformation u = 1/v equation 
(25) becomes 



2.d'P „ dP 2,d^P 



(26) 



We use the technique of separation of variables; let 
P{u,iJ.) = R{u)M{fi) and dash denotes differentiation with 
respect to u, equation (26) becomes 



(1 - u^)R" - 2uR' _ (1 - 



-z(/ + i) 



(27) 



R M 

Equation (27) for R(u) is the Legendre differential 
equation, therefore the solution is a linear combination of 
R{u) — CaPiiu) + CbQiiu), that has to be finite in the inter- 
val u G [1,00), therefore C5 = 0, as Qi(u) becomes infinite 
at u = 1. 

The angular part ca n be solved analytical ly using the 
following transformation (|Wolfson fc Lowlll992l ). M = (1 — 
/i^)^^^, therefore the second equality of (27) becomes: 



(1 - At^)A^r " 2/iA/( + l{l -f l)2Afi = 0. 



(28) 



This is again the Legendre differential equation, and the 
solution for the angular part of P is M{jj,) = (1 — fj?)Pl{^). 
Therefore the ffux function of the field is 



>,m) = { 



Fmax{l - iJ^)Pi{lJi)Pl{u) U ^ 1 







u <l 



(29) 



3.1.1 The vacuum expanding dipole 

Let us consider 1 = 1, that corresponds to a field where 
a dipole of magnetic moment M = Mpti is i n itiate d at 
t = Q. Following the formalism of IS ommerfeldl (Il952l ) for 
the Hertzian dipole modified for a magnetic dipole we can 
express the electric and the magnetic fields using the quan- 
tity II which is a vector pointing at the axis of dipole z and 
is a function of position and time; dot denotes derivation 
with respect to ct 
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E = -V X n, 



(30) 



spherical shell c{t — to) < r < ct and (III) r > c(t + 
D{ct — r) becomes: 



B = -n + v(v-n) (31) 

It is shown in appendix A that for any function D{ct—r) 



(32) 



and the fields satisfy Maxwell's equations. The fields 



D ■■ 



m{ct — r) 




7 + 



t < 



ID. 

2 

IQ. 

2 

r TQ^ 

c 2 



< t 

< t < 



+ 



(I) 

ill) 

{III) 



(39) 



This substitution iucrcascs the complexity of the formu- 
lae for the fields but removes the singularity at the horizon 
and allows us to study in greater detail the structure of the 
field near the horizon and how it depends on the time-scale 
for the dipole to establish a linear behaviour. The fields cor- 
responding to D are given by putting D into (33) and (34) 



E : 



sine / " D\ 2 



(33) 



2 cos I 



\^^-y^—(p^-^^Y (34) 

Therefore, we can evaluate the electric and the magnetic 
fields for a linearly increasing dipole that is switched on at 
f = and lies at r = 0, it corresponds to D = mict — 
r)H{ct — r), where H is the Heaviside function which is 
unity for argument greater than zero and zero otherwise. 
We divide the space into three regions. The first one (I) is 
r < ct; the second (II) is the surface r = ct and the third (III) 
is r > ct. In the first region there is an electromagnetic field 
because of the dipole; in the the third region the field is zero, 
as the message of switching on the dipole has not arrived 
yet; the second region is the horizon surface of discontinuity. 
Therefore in region (I) the fields are: 



Br 



Sine ( 

CTQV \ 



2 + 



2(e(t+^)-r) 



2mct^ 

2m cose / {c(t+^))^ 



,0 ( { 



m sin 







{I) 

in) 

(III) 



{!) 

-r) (77) 
(777) 

(I) 
(777). 



(40) 



(41) 



(42) 



These are perfectly regular and acceptable. 
The higher order self-similar multipoles are given in Ap- 
pendix B. 



E : 



sin^ 



ct>, 



B = mc^ {2 cos er + sin 90). 



(35) 



(36) 



Since the fields are related to derivatives of D and the 
first derivative of D is discontinuous at r = t we expect 
singularities on the surface r — ct (II). Indeed there are sin- 
gular surface fields. These infinities we should have expected 
because in the analogous electric dipole case velocities are 
suddenly imposed on the charges to make the linearly grow- 
ing dipole. Thus initial accelerations are infinite so the power 
radiated should be singular. 



E = -m 



sin 6 



(^Sict - r) 



+ 



H(ct 



(37) 



mc / 2t cos 9H{ct — r) ^ 



r + 



sme(S{ct -r) + 



tH{ct - r) 



(38) 



These fields are unacceptable. In region (III) there are 
no fields at all. To avoid the singularities in (II) we take into 
account the finite time needed to accelerate the dipole to a 
given growth rate. To do this we smooth by averaging the 
fields over a small spread to of switch-on times. Thus region 
(I) now extends r < c(t + region (II) spreads into a 



3.2 Solutions containing currents 

In the previous sections we have given analytical expres- 
sions consisting of magnetic fields that lie on the meridional 
planes and electric fields that axe azimuthal; such fields do 
not contain currents, or charge densities and each arises from 
a growing multipole at the origin. In this section we seek 
analytical solutions for electromagnetic fields in which the 
magnetic field has an azimuthal component, the electric has 
a meridional component and there are charge and current 
densities in the space. 

We need to solve (10), Prendergast's equation. We seek 
a solution that is a product of a function depending on v and 
one depending on n; using the usual transformation v = 1/u 
equation (10) reduces to 



d_ 
du 



We multiply by {v? — 1) /P and look for solutions in the 
form of P = g(u)F[ij) 



g du 



I 2 ,.dg 



F diJ? 



(44) 

PdP ^ ' 



We set: 
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U{u) = 



1 d 



du 



- d^F 
F dy? ' 



PdP 

so the equation takes the form: 



U{u) + {u^ - l)M{p) = -Q{P), 

operating with F/F' d/djj. we find 



(45) 
(46) 
(47) 

(48) 

(49) 

C{gF), with 

A{g{u)) = u'^-l and C{P) = -PQ'. Taking logs and oper- 
ating with gd/dg we find that (i) d In j4/dln <; = din C/dIn P 
or alternatively (ii) C = B = which leads to Prendergast's 
linear model P{P) oc P. Henceforward we consider only case 
(i), F is not constant since P can not be a function of g 
alone, so both d\nA/d\ng and d In C/d In P must equal the 
same constant 2/1. This is a generalisation of the I used be- 
fore in the current free case. Thus without loss of generality 
g = ^'/2 ^ (u^-iy/-^ and C = -PQ' = ~KP^^' where K is 
constant. Finally from (77) FM'/F' = -KP^'''/{u^ - 1 
-KF^^\ hence: M' 
Q' = KP'^'^-^ so Q 
now reads: 



t 2 ,.FM' 



which is of the form A{g)B{F) 



-KF^/'-^F', M = -i^KlF^^' + c), 
= iif/P^/' -f c', our basic equation (44) 



/ 2 ixl-i/2 d 



/ 2 ,x " / 2 ,Ni/2 



-{u^ - 1){\kIF^" + c)= ~{\kIP^" + c'), 
but 



(u2_i)i-i/2Ar 

du 



Z(m^ - l) + ru 



(50) 



(51) 



so equation (50) reduces to a linear equation in u^. Com- 
paring coefficients of 



1(1 + 1) _ {^KIF'''' + c) = -]^KIF^'' 

and from the coefficients of vP in (50) 



Z+ ^KIF^^'^ + c= ^KIF'^^^ 



(52) 



(53) 



Evidently we must choose c — 1(1 + 1), so the first is 
satisfied and c' — —l^ so the second is too. However in the 
non-relativistic approximation v? » 1 so then the second 
equation is all negligible compared with the first. 

From the definition of M we have: 



(54) 



(1 - f/)F" + 1(1 + 1)F = ]^KlF^+^'\ 

and by normalising / = F/Fmax, so / ^ 1 

1 



(1-m')/" + Z(Z + 1)/ 



2/1 rl + 2/l 



(55) 



(56) 



which is the same equation as that studied by 
iLvnden-Bell fc Boilvl l|l994^ ■ who found it in non-relativistic 
MHD. We choose the solutions obeying boundary conditions 
that / is zero at 9 = 0, [fi = 1) and at any chosen outer 
boundary Q, (/io). The constant C = ^KlF^a^ is then de- 
termined in such a way that / is one at its maximum. For 
small I and — 1 — ^6^ the solution for / is 



1 20^ 
/=1_ Incosh 



where v = 1 + l/l and g is close to 1 and given by 



(57) 



1-VT 



-2v , 



g — — cosh ^ — \ -\ — In [l 



~ H- -ln2. 

The value of C associated with this solution is 



(58) 



(59) 



i^;^2/! ^ ^ ^ 
2 //92 

Equation (57) gives solutions which are symmetric in 9^ 
around their maximum e^/2. This is the case for I = Q, when 
I is small but not zero the maximum is slightly displaced. A 
more accurate formula taking into account this displacement 



/ = 1 - i Incosh (4^{fJ. - m 
V \ sm 9i 

m cosh|^ — —{^. — MijJI 



sm V\ 



(60) 



75 6(i^+2in2) ' A*! — 1 — \9 \ and C = C^v — l^v, for details 



where u — v 

e 1 



on the derivation see iLvnden-Belil (|2006l ) (appendix). 

The remaining equation comes from the definition of Q 

(47), 



2{l + l) 



32(1-1-1/0 



ZVVC2, 



(61) 



where C2 is constant. If Z > then has a minimum at 
Po = {2l/Ky^^. We choose C2 so that this minimum is zero 
and take /3 to be zero for P < Po so 



13' 



I /2Po2{^[(P/Po)^+^/'-l] - [(P/Po)^-l]} 
\ 0,P<Po 



(62) 



Had we chosen a lower value of C2 then (3'^ would have 
been negative for values of P near Pq. However must be 
positive. An unfortunate consequence of choosing the upper 
expression (62) for for all P is that / when P = 0. 
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This entails an infinite near the axis which corresponds 
to a hne current there. To avoid this in a continuous and dif- 
ferentiable way we choose the upper form (62) only where 
P > Po- Elsewhere we choose (3 = 0, the current-free so- 
lution. Wo avoid the use of the linear form /? oc P because 
those solutions have an infinity of oscillations in the sign of 
P as they approach u = 1 (see Prendergast 2005). So they 
contain fields that are disconnected from the origin. 

Thus, our solution of the relativistic Prendergast equa- 
tion is 



P = P™a./(M)K-l 



,V2. 



P^ Po 



(63) 



where / is given by (57) or more accurately by (60). 

In what follows we determine the fringing field in the 
cocoon region P < Pq. We choose the boundary to be a cone 
with a total opening angle (not semi-angle) of less than 22°. 
This is the case of greatest importance for applications to 7- 
ray bursts, however we should point out that exactly similar 
methods work for any cone or even for the whole sphere. It is 
just that Legendre polynomials or more generally Legendre 
functions of replace the simpler Bessel functions of the 
narrow cone cases. 



3.2.1 Solution in narrow cones 

Since our main applications will be the jets and 7-ray bursts 
it makes sense to study solutions within conical boundaries 
with quite narrow cones 9 ^ Q. Provided the total opening 
angle (not the semi-angle) of the cone is less than 22° the 
actual opening angle does not affect the form of the solu- 
tion whereas at larger angles it does. Thus there is a real 
advantage in studying the small angle case which allows us 
to replace high order Legendre functions of angles by Bessel 
functions of order 1. The 11° limit comes from our replace- 



ment of (1 + /i) = (1 + cos f 



2- 



by 2. The fractional 



error in neglecting ^6^ as compared with 2 is which is 
less than 1% for 6* < 1/5 ~ 11° 

In the region where P < Pq, 0^ is zero so the solutions 
of (43) that we need, have no currents, however they are 
solutions with a very awkward boundary on P = Pq. In 
place of trying to fit a current-free solution on that boundary 
to our known solution within the region with P < Po we 
adopt the current-free solutions within the whole cone as 
our complete set. We expand the known right hand side of 
equation (43) in terms of our Bessel functions of angle within 
the cone. We then get a known function of w as a source 
term on the right hand side of each Bessel component of the 
equation (43). We solve for the radial functions of u which 
are the coefficients of the different Bessel components, each 
with its own known source term. 

Under the assumption of narrow opening angle and in 
the absence of currents the angular part of (43) becomes 



setting J = 9^ the above equation becomes: 

■1 d /„dJ^ 



9dr9{'d9)+''^'' + '^'. 



0, 



(64) 



(65) 



-9T9['-9)+''^' + ^^'-'^''- 



(66) 



So far, we have only defined J up to an arbitrary ad- 
ditive constant which is eliminated in J so we absorb the 
constant into J and get: 



1 d ( dj 



(.^)+K^ + i)J^o. 



(67) 



The solution is aJo{k9) +bYo{k9), where + = k^. 
If we want a solution regular at = we must omit the b 
term so J = aJo{k9), and 



j = ak9^M^ = -azMz), 



(68) 



d{k9) 

where z = k9. 

In order that J should vanish at O we need Ji(fcO) = 0. 
The solutions for k are ks = js/Q where ji = 3.83, j2 = 7.01, 
js = 10.17 etc p 409 of A & S, and je is the s-th zero of Ji {z) 
for z > 0, called jis in A & S. 

Now any function F{9), F(0) = F(e) = can be ex- 
panded as a series. 

F(e) = ^F.Ji(i.e/e). (69) 

s 

The orthogonality relation is A & S 11.4.5 



tJl{jst).Jl{j„t)dt = - [J[{jn)]'^5s 



(70) 



To find Fs we multiply by 9/QJi{jn9/Q) and integrate: 

[|ji(jne/e)^F.Ji(i.e/e)]d(|) = 



(71) 



So the F„ corresponding to any F can by found by 
integration and a knowledge of the j„ and Ji{jn) found in 
the A & S p 409. 

Now write: 



and 



(72) 



(73) 



Then multiply (43) by Ji{j„9/Q) and integrate from 
to 0; on division by ^[J[{jn)]^ we get, since {9/@)Ji{js9/Q) 
satisfies (67); 



d_ 
du 



du 



||p„(M)=ir„(«). 



(74) 



The above equation is Legendre's with a known right 
hand side K„{u). This may be solved using the boundary 
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conditions by variation of parameters. We write each p„(u) 
in the form 



(75) 



where v{n){v{n)-\-l) — j'n/Q^ and P^, Qi, are Legendre 
functions satisfying 



d 

du 



[l-u)-r^ +v{u + l)P^ = Q, (76) 
au J 

and Qn likewise. We choose An and Bn to satisfy 



X(w)-P^(n)('i) + B'^{u)Q^(^n) = 0. 
Thus 

p'„ = A„PI{u) + B„Q'^{u). 



(77) 



(78) 



Substituting (78) into (74) the terms with A and B 
undifferentiated vanish since P,j and Q„ satisfy Legendre's 
equation so the only terms that survive are 




Figure 1. The fiux function P for u = 1.5 and I = 0.1, given 
by equation (63) and reconstructed by using equation the first 
15 components of (100). The small difference of the two forms is 
plotted at the bottom. The deviation of the two docs not exceed 
0.2% of the maximum. The part of the plot between the vertical 
dotted lines corresponds to P > Pq. 



(1 - u^MnfPl + BnQ'A = K„iu). 



(79) 



We solve (77) and (79) using the Wronskian W = PQ' - 
Q'P. But the Wronskian is -(u^ - 1)"^ see A & S 8.1.9., so 



(80) 



Aniu) — A„(ub) = / KnQu(n)du 



similarly, 



Bn 



KnP„(n)du 



(81) 



where Ub is the value of u at the surface where the 
field expands at some given time, and ut is the value of 
It at highest point of the cocoon P + Pq. An{ub) can be 
determined by analysing P{ub,6) into Ji and by (75) 



An{ub)Pv{ub) + Bn{ub)Qu{ub) = Pn{ub) 



(82) 



One must choose the end points of the integrations to 
ensure the boundary conditions are satisfied. Indeed formula 
(72) provides accurate results when the first few terms are 
evaluated and we are not very close to the extreme relativis- 
tic limit. In the paradigm solved for Z = 0.1 the deviation 
between the reconstructed solution and the solution satis- 
fying (56), which holds inside the lobe where P > Po, is 
0.2% at u = 1.5 and 0.3% at u = 1.05, 0.8% at m = 1.01 
and 4% at m = 1.001 by taking into account the first fifteen 
terms (Figures 1-3). As u comes closer to unity the deviation 
between the reconstructed solution and the self-similar one 
increases. This deviation is systematic as the reconstructed 
solution always gives a flux function that is greater than the 
self-similar one. This is because the self-similar solution is 
proportional to — 1 and becomes at u = 1, whereas 
the Legendre P functions converge to 1 at the top. Very 
close to the top of the lobe we reconstruct the solution by 
analysing P{ut,9) into Bessel functions, then by virtue of 




Figure 2. The flux function P for u = 1.01 and I = 0.1, given 
by equation (63) and reconstructed by using equation the first 15 
components of (100). The difference of the two forms is plotted 
at the bottom. The deviation of the two does not exceed 0.8% of 
the maximum. The part of the plot between the vertical dotted 
lines corresponds to P > Pq. 



(76) and setting Bn{ut) — because of (81) we can deter- 
mine An{ut). Then from (82) it is clear that for 1 < u < ut 
An{u) — An{ut). Thus we have a complete description of 
the field everywhere. 



4 PHYSICAL QUANTITIES 

In this section we evaluate physical quantities such as the 
twist, the pressure and the energy corresponding to the 
mathematical solutions we describe above. We use quantities 
appropriate for 7-ray bursts of the long-soft type. We study 
two cases; one with semi-opening angle of 0.2rad which is 
the maximum value for the assumption of narrow cones to 



Fields from a relativistic magnetic explosion 9 




25 



20 1* 



15 



E 

© 10 



10 

I 



15 



20 



Figure 3. The flux function P for u = 1.005 and I = 0.1, given 
by equation (63) and reconstructed by using equation the first 15 
components of (100). The difference of the two forms is plotted 
at the bottom. The deviation of the two docs not exceed 1.3% of 
the maximum. The part of the plot between the vertical dotted 
lines corresponds to P > Pq. 



hold and one with semi-opening angle of O.OSrad which is 
closer to what is expected for 7-ray bursts. We find that 
in the latter case the fields are more twisted so the energy 
of the initial purely poloidal field is amplified by a greater 
factor. 



4.1 Twist 

The form of the field defines the distribution of the twist. 
For the reasons discussed in the section 3.2.1 we use the the 
solution for the fiux function given by (63) for P > Pq. We 
have taken as a given that the magnetic field is contained 
within a cone. Although a detailed hydrodynamic study is 
needed in order to find the shape of the magnetic cavity; 
a pressure environment that depends in r like p oc r~'^'~'* 
leads to the conical jets. 

Inspection of (62) shows that toroidal torque f3 increases 
with P when P > Pq but it is zero for P ^ Pq. However 
our solution (63) shows that the lines of force for large P 
do not go to large radii where the twist is. Thus the twist 
of the lines is greatest for some P — Pm greater than Pq 
and less than Fmax- When I is large Pm is close to zero, 
but as I decreases it shifts closer to Fmax- The maximum 
twist increases as I decreases (Figures 4 and 5). Smaller I 
corresponds to fields that decrease slower with r and at the 
same time they demand greater twist. 

The lines of force of the magnetic field are given: 



dr 

Br 



rdO r sin 
Ba B, 



4> 



(83) 



In the region where P < Po there is no twist in the 
field, as f3 — 0, thus B^ — 0. In the region where P > Pq 
the field can be evaluated by virtue of (63) and (8) and by 
using the second equality of (83) the total twist $f for a 
field line determined by P is given by: 



Figure 4. The maximum twist with I for a cone of semi-opening 
angle of 0.2 rad. As / decreases the maximum twist inside the jet 
becomes greater. It becomes zero at I = 18.6. 
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Figure 5. The maximum twist with I for a cone of semi-opening 
angle of 0.03 rad. As I decreases the maximum twist inside the 
jet becomes greater. Compared to figure 4 the twist is bigger, it 
drops to zero for I = 127.2 which is not included in this plot. In 
general the qualitative behaviour of #mai(0 in similar to this of 
the wider cone. 



$P = 



IP 



d9 



(84) 



where 9i and 62 are the polar angles of the points where 
the field line rises and sinks at the base. We then express u 
as a function of 6 for some given P and we integrate. This 
gives the detailed shape of the field lines (Figure 7). Most 
of the twist is concentrated near the top of each field line. 

From integral (84) for 9i and 62 corresponding to the 
footpoints of the field lines we evaluate the total twist a field 
line carries. The maximum twist depends on the opening an- 
gle of the jet. The overall behaviour of ^max can be approx- 
imated by the analytical formula for the total twist in the 
non-relativistic case (jCourgouliatosI [ioOSh 4^ (Figure 7) 
or -^^^ for the case of narrow cones we are interested in this 
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Figure 7. The maximum twist as a function of the semi-opening 
angle for I —* 0. The maximum twist is inversely proportional to 
the semi-opening angle, but instead of proportionality constant 

y/2lT = 4.40 it is ^max 



4.03 
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Figure 6. Left: Field lines in the case of ( = 0.1 in the jet of 
semi-opening angle 0.2rad. The field lines correspond to fluxes 
0.6 red and 0.7 blue of the total flux. The outer pink cone (light 
coloured in black and white) is the cone confining the jet, whereas 
the blue surface (dark coloured in black and white) is the surface 
corresponding to flux 0.6 on which the blue fleld line lies. 
Right: The lobe and the cocoon for the same jet. The lobe in blue 
comes very close to the top of the jet. The size of the cone in both 
cases is ij» . Narrower jets have many more twists. 



study. Thus the smaller the opening angle the more twisted 
the field is, giving evidence for the relation between colli- 
mation and the number of the turns performed. Although 
the total twist is inversely proportional to the semi-opening 
angle, the constant of proportionality is somewhat smaller 
than what was found from the non-relativistic analytical for- 
mula. 

The total twist for a given field line does not de- 
pend on time, indeed they merely expand with the ex- 
panding frame. Nevertheless the field lines rotate as they 
expand relative to the stationary frame. The azimuthal 
velocity component given in equation (17) describes that 
rotation of the field lines and is a fraction of the speed 
of light, its maximum value is 0.6c. As we have a non- 
stationary solution of Maxwell's equations, problems at 
the light cylinder cannot occur, and even a radius for 
the light cylinder is not defined. Rotation has been ob- 
served in non-relativistic jets emerging from T-T auri stars 
llCoffev. Bacciotti. Rav. Eisloffel. fc Woitasll2007[ ). 



4.2 Energy 

We can easily find the pressure and the energy carried by a 
given electromagnetic field since we know the actual fields. 
As the energy density of an electromagnetic field is propor- 
tional to the pressure we are going to treat them together. 
The energy density is given by 



e=i-(B^ + E^). (85) 



4-2.1 Energy distribution 

From equations (8) and (9) and the flux function P from 
equation (63) we can evaluate the pressure at some given 
time t (Figures 8 and 9). The energy is mainly concentrated 
near the base of the jet and near the top of the jet. It is 
expected that the energy is strong near the base as the fields 
decrease with distance from the centre. However unlike other 
cases studied where the fields extend to infinity, there is a 
very strong field near the top. This is because of the turn 
over of the field lines, which are confined inside the light 
sphere, so they have to close. Thus, they have a very strong 
Bg component. 

As we approach the light horizon where the expansion 
velocity is equal to the speed of light the energy strongly 
depends on the smoothing of the fields. If we considered 
an infinitely fast turning on of the magnetic field then it 
would carry infinite energy, however this is unphysical as it 
takes finite time to switch on. That time scale determines 
the energy of the initial pulse; it is inversely proportional 
to the time-scale of switching on the dipole. The physical 
system we have in mind is the following: as the core of a 
massive star collapses to form a compact object or a black 
hole it spins much faster and winds up the magnetic field 
where it connects to the rest of the star or to any disc that 
may form. This leads to a dramatic increase in the energy 
by two orders of magnitude in the case of a narrow jet, 
as shown in the next section. All the twist is imposed in 
this stage. The concentrated energy of the field leads to an 
explosion. The magnetic field expands in the star and a jet 
forms. The scales chosen for this study are J?, — 10^'' cm for 
the size of the jet. This is the typical size of the jet related 
to the prompt emission of a 7-ray burst. Rc = lO'^cm for the 
smoothing scale that is a few times the length scale of the 
Schwarzschild radius of a black hole of IOMq . 
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Figure 8. The energy density and the static pressure in a section 
of the jet. The contour lines correspond to surfaces of constant 
pressure and differ with the consecutive by a factor of two. They 
are stronger at the bottom and they have a hot spot near the 
top. This is because the field lines turn back near the top giving 
a strong Bg component in the field. They increase as we move 
towards the centre of the jet reaching a maximum and then de- 
crease in the outer part because of the strong field inside the 
lobe. The pressure near the axis is stronger than the pressure near 
the edge because of hoop stresses. The axes are not to scale, the 
units of the axes are fraction of the total jet length at that time. 



Figure 9. The energy density and the pressure in 9^ and r. The 
same conclusions as from figure 8, but in this case one can see the 
symmetry of the structure in 6^. 
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We consider a magnetic field of order _Bo = 10 — 
10^*Gauss at Ro = 10^ cm powering the jet. This magnetic 
field comes from the central compact object, a discussion on 
the origi n of strong magnetic fields in neutron stars can be 
found in lSpruitI ()2008 ) and references therein. The strongest 
magnetic fields estimated for magnetars are lO^^Gauss. We 
assume that such a field originates at Rc and then expands 
outwards, as the surrounding pressure near the centre drops 
relatively slowly. A constant pressure environment leads to a 
field which is confined within a cylinder without decreasing 
at all, we take an intermediate case for which the field be- 
tween Rc and Ro drops by two orders of magnitude. Above 
this radius and below R, we apply the self-similar solutions 
inside a narrow cone. These give a total energy of order 
W^^erg in the jet from which 10*^ erg is due to the switch- 
ing on of the field concentrated in a very thin layer just 
below the fight sphere. The width of this layer is lO^cm and 
the fields contained are of order lO^Gauss whereas the field 
at the main body of the jet at 5 x 10^^ cm is of order of 
10^ Gauss. The energies quoted above and on the tables are 
at t — 300s after the initial explosion in the centre of the 
star, when the top of the jet reaches the surface of the star. 



Table 1. The maximum twist and the energy of the various 
components of the magnetic field for a given flux Fmax = 
2 X lO^^Gauss x cm?, corresponding to an average magnetic field 
of order lO'^^ Gauss emerging from a spherical cup of opening 
semi-angle 0.2rad and radius 10^ cm. As / decreases the field gets 
twisted and its energy amplified. 

4-2.2 Energy as a function of twist 

The total energy depends on / (Figures 10 and 11). It can 
be seen (tables 1 and 2; figures 4 and 5) that I parametrises 
the twist. We prefer the use of the twist as parameter as it 
describes a physical process rather than I which is just an 
index. Thus a decrease in I can be viewed as an increase on 
the difi'erential rotation at the base. When difi'erential rota- 
tion is imposed at the base there is a toroidal component 
in the field which increases the energy of the configuration, 
and in addition to that the poloidal part of the field ex- 
pands to reach a force free state. This dual process leads 
to an increase in the total energy. The factor by which the 
total energy is amplified depends on the maximum twist 
the configuration can tolerate and thus the opening angle 
which determines the maximum twist. Narrower cones get 
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Table 2. The maximum twist and the energy of the various 
components of the magnetic field for a given flux Fmax = 
10^^ Gauss X cm^ , corresponding to an average magnetic field of 
order 10^* Gauss emerging from a spherical cup of opening semi- 
angle 0.03rad and radius lO^cm. As I decreases the field gets 
twisted and its energy amplified. This case clearly illustrates that 
a jet of smaller opening angle can be more twisted than one of a 
bigger opening angle. The I for which the field does not have any 
toroidal component is 127.2. The energy of the purely poloidal 
structure is much smaller than this of the one that is twisted. 
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Figure 11. The energy as a function of the maximum twist for a 
cone of semi opening angle of 0.03rad. The conclusion are similar 
to those of figure 10. The qualitative behaviour of energy with 
twist is the same for smaller opening angles. 
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Figure 10. The energy as a function of the twist for a cone of 
semi opening angle of 0.2rad. It is clear that as the twist increases 
energy is injected in the system. Initially all the the energy is 
in the poloidal field, however the energy of the toroidal field, 
increases reaches a maximum and then decreases. Unlike other 
cases of spherical symmetry where the field lines are allowed to 
reach infinity there is some contribution of the toroidal magnetic 
field to the total energy of the field when it reaches its maximum 
twist. 



a greater amplification on their total energy; in the limiting 
case which is equivalent to an infinite cylinder, the twist has 
no constraint, thus the energy can be infinitely amplified. 



4.2.3 Energy flow 

As the whole structure expands there is a flow of energy in 
the jet. The energy contained in an element of the the jet 
from angle 9 to 9 + d9 and in velocity space between v and 
V + dv is 



El 



1 



rdP 



IGtTT^u^ sin 9 



This energy is inversely proportional to time. The 
Poynting vector for the fields is 



S = 



cE X B 

47r ' 



(87) 



substituting from (8) and (9) 



S = 



fl- 
op dP- 

U 



(dP_y 



\ dv J 
vT dP - 



The 9 component of the Poynting vector is zero at the 
boundaries of the jet and on the axis, so it is consistent 
with the demand that the jet is constrained within a cone 
and there is no fiow of energy out of the sides of the cone. We 
now apply the Poynting theorem for the expanding volume 



8n 



jL 

cdt 
v-dS 



E' 



8-K 
E X B 

47r 



dV = 



dS. 



(89) 



As opposed to the usual form for the Poynting theorem, 
there is an extra term in the right hand side because of the 
expansion of the volume. The term in the left hand side is 
the integral of equation (86). The terms in the right hand 
side show in which direction the energy is transfered in v 
space. The integrals should be performed on the surfaces of 
the inner and the outer boundary of the expanding volume. 
Nevertheless, by integrating only on one surface we can find 
whether the energy flow is towards larger or smaller v. 



E X B 



(8^) 



47r 

[ ^ 2(dP.2.1~V^ l,dP.2\d9 



dS- 



E^ + B' 



V- dS = 
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The integral of equation (90) is negative when inte- 
grated for the solutions of P and T found in section 3.2. 
Thus the energy flows towards smaller v. Nevertheless the 
energy flow relative to flxed axes is outwards everywhere. 
This discrepancy occurs because the energy flow as described 
by the Poynting vector is slower than the velocity due to the 
uniform expansion. 

The region at the top has to be treated separately. Near 
the top r = and it is only the component due to ^ 
that contributes to the energy. The energy in the top is 
compressed in a thin shell (region II). The width of this shell 
is Re- Thus near the top the flux function drops from P{v = 
1,6) to within a distance Rc as described in Appendix 
B for the smoothing of the flelds in the case of self-similar 
multipoles. P{v — 1, 9) is approximated with good accuracy 
by Pof{9), where Po is defined for equation (62) and f{9) is 
the solution to equation (56), note that the exact value of 
P{v = 1, 9) is not given by equation (63) but it is a sum of 
multiples of Bessel and Legendre functions. After we smooth 
the field we find that the energy contained in region II 



Eu 



1 Pi 



Atv^ Rc /(, smf 



d9. 



(91) 



Thus the energy of the shell (region II) is constant. 



4.3 Pressure 

4-3.1 Static pressure 

As the pressure of the magnetic field is proportional to 
^-4-21 £^j. £j^g ^^gg q£ ^j^g self-similar solutions in the non- 
relativistic region the field shall demand a special pressure 
environment so that it is contained within a cone. Let us 
assume that the pressure of the surrounding medium is de- 
scribed by a power-law r"'', thus as the field is differentially 
rotated it will expand primarily up the axis and if the in- 
dex of the pressure of the medium <^ is greater that 4 — 2/ 
then the jet will not have a constant opening angle but will 
expand sideways. 



4-3.2 Stagnation pressure 

Near the top the case is rather different, as there is a partic- 
ularly strong ram pressure. This is because the head of the 
jet encounters the gas of the medium in a velocity Vh that 
is close to the speed of light, let this velocity correspond 
to some Lorentz factor 7. As there is only a Bg component 
there, in the comoving frame the magnetic field is: 



B'g 



(92) 



This magnetic field encounters an environment of den- 



sity 



P =1P, 



(93) 



where p is the density of the medium in its rest frame. 
The 7 factor is because the coordinate in the direction of 
expansion is shrunk by a factor of 7. Balance of momenta 
gives a stagnation pressure in the comoving frame 



'2 



/ 2 2 B(, 



= 7 



8n' 



(94) 



Setting the stagnation pressure equal to the magnetic 
pressure of the top field we can determine the 7 that the 
magnetic field reaches as it expands. This depends on the 
magnetic field and the density the jet encounters; as Bg and 
p are known we can easily find the value of 7. For densities 
around 10~^ g/cm^ a Lorentz factor of 2 is found,whereas for 
densities smaller than 10~^g/cm^ Lorentz factors greater 
than 10 are found, and for the case of an HII region sur- 
rounding a massive star the jet becomes extremely relativis- 
tic with 7 >> 100. However, near the centre of the star the 
density greatly exceeds the value of 10~^ gr /cm'^ , thus the 
jet pushes its way through the material at non-relativistic 
speeds. However it is possible that a funnel forms by an other 
mechanism. In the case of 7-ray bursts there is a black hole 
formed in the centre of a rotating star thus, the material on 
the axis will fall into the black-hole. Then there will be an 
underdense region o n the axis allowing the expansion of the 
jet, see e.g. iLvnden-Bell (,197a l. 



5 APPLICATIONS 

Relativistic jets were first found emanating from radio galax- 
ies and quasars, then more spectacularly still from 7-ray 
burst sources and microquasars. Related but non-relativistic 
jets arise from young stars with accretion discs and generate 
Herbig-Haro objects. We consider that all these phenomena 
come from the magnetic energy generated when a magnetic 
field is wound up by the differential rotation of its foot- 
points. The opening angles of many of these jets is of order 
4°, giving semi angles of 2° or so. 

It is believed that massive Wolf-Rayet stars in rapid 
rotation generate the long-soft 7-ray bursts when their cores 
collapse event ually to make black holes and the associated 
supernovae (|Heger. Fryer. Wooslev. Langer fc HartmannI 
,2003 ). Then the differential rotation is between the col- 
lapsing core of the star and the enclosing envelope. It has 
also been suggested that the outco me of the collapse is a 
strongly magnetised compact object (|Duncan fc ThompsonI 
Il992h releasing energy in magnetic field. This magnetic 
energy is generated by differential rotation deep within the 
star. While the strong magnetic field expands to balance 
the external gas pressure the density of the external 
envelope initially prevents relativistic expansion along the 
rotation axis. Under expansion at constant flux the field 
in a magnetic bubble or cavity varies, with length scale, 
L, as and the magnetic pressure as L~*. Thus while 
it is within the body of the star the magnetic cavity is 
coUimated when the ambient pressure falls less rapidly 
than r~* subtending smaller angles at the star's centre 
as r increases. Such configurations of non-re l ativist ic jets 
were described in ISherwin fc Lvnden-B"eiil and 
references therein and ha ve bee n applied to 7-ray bursts by 
Uzdenskv fc MacF adveni l|2006t ). In the outer parts of the 
star the pressure falls more rapidly than and the head 
of the magnetic cavity accelerates until it breaks through 
the stellar surface and expands relativistically. Numerical 
studies of ultrarelativistic MHD jets have been done by 
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iFendt fc Quvedl { 20041 ) ■ where it was found that high 
colhmation angles can be achieved by a configuration where 
the toroidal component of the magnetic field dominates. 
When the central object is a rotating bl ack hole then the 
Kerr geometry was taken into account (|Fendtlll997l '). In 
the case we study the main source of confinement is ram 
pressure on the low density stellar wind. 



Jets of narrow opening angles can be more twisted than 
those of wide opening angles, giving evidence that colhma- 
tion increases with twist. 

Finally a strong stagnation pressure is found at the head 
of the jet. It prevents the jet from being relativistic near the 
base but as it expands outwards and encounters underdense 
material it becomes extremely relativistic as is indeed ob- 
served in 7-ray bursts. 



6 CONCLUSIONS 

In this paper we have solved analytically the equation pro- 
posed by Prendergast for relativistically expanding axisym- 
metric self-similar force-free fields. This equation is a rela- 
tivistically moving extension of the Grad-Shafranov equa- 
tion for static force-free fields with axial symmetry. 

To excellent accuracy our solution is given everywhere 
by equations (62) and (63) with the /(^) function given by 
equation (57) or more accurately (60). The resulting elec- 
tromagnetic fields are given by equations (8) and (9). All 
solutions of the Prendergast equation need to be averaged 
over a small time At that corresponds to the light crossing 
time across the magnetic explosion. His equation assumes 
a point explosion. Exact solutions of Prendergast's equa- 
tion give singular fields where the expansion speed is that 
of light but we have shown how averaging over Ar yields 
genuine solutions of Maxwell's equations without singular- 
ities. Nevertheless this Ar may be quite short and the re- 
sulting fields close to r = ct are strong, carrying around 1% 
of the system's energy. This may be related to 7-ray burst 
precursors. However, the fields at the top require further 
modification as the ram pressure of their relativistic mo- 
tion is retar ded by the interstellar m edium. Considerations 
as those of (|Rees fc Meszaro^ |2005| ) are clearly important 
there, whatever drives the shock waves and it is there that 
the eponymous 7-rays themselves are generated. 

While this paper has given a basic mechanism that pro- 
duces highly relativistic motion there is a big gap between 
that and the phenomenology of 7-ray bursts. 

On the way to understanding the solutions of the fully 
relativistic Prendergast equation we were led to give purely 
poloidal solutions of Maxwell's equations for time-dependent 
multipoles, interesting in their own right. We have also 
seen how in the non-relativistic limit our pr o blem is re- 
lated to the solutions found in iLvnden-Belll (|2006D and 
iLvnden-Bell fc Boilvl l| 19941) . 

Finally we have seen how small opening angles of the 
cone correspond to greater coiling of the magnetic field and 
thus to mush greater concentration of energy into the jet. 

Setting some values for the magnetic field quoted from 
the stronger magnetars we find that we can power jets up to 
the energies of 7-ray bursts, for which the isotropic energy is 
between 10^^ to 10^'' erg (Piran 2005). In the case studied we 
found that the energy inside the jet is of order 10^^ erg, where 
the isotropic energy is two orders of magnitude greater. The 
very strong magnetic field found at the top of the magnetic 
configuration can oppose the strong ram pressure because of 
expansion and swipe material out of the way for the jet to ex- 
pand. This feature may b e related to the precursor observed 
in som e 7-ray bursts e.g. iMurakami.Inoue. Nishimura et aU 
(|l99ll ). as it is the first to break out and carries a small 
fraction of the energy of the 7-ray burst. 
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APPENDIX A: DIPOLE 

Following the formalism of ISommerfeldl (|l952h we express 
the electric and the magnetic fields in the following way, dot 
denotes differentiation with respect to time, 

E = -V X n, (Al) 



B = -n + v(v-n) 



(A2) 



The fields must satisfy Maxwell's equations in vacuum. 
By construction they satisfy 



V-E = 0, 



(A3) 



n 



(AlO) 



where V • A = and A = 0. Using the gauge transfor- 
mation that leaves the fields unaffected as we have shown 
above, we can substitute in favour of II, and equation (AlO) 
becomes: 



n-v^n = A-v^(Vx c) 



(All) 



Note that the divergence of the RHS vanishes. Thus by 
choosing 



V X C = 

4-K 



(A12) 



we can make the RHS of (All) vanish and therefore, 
equations (A5) and (A6) reduce to the wave equation. 

The next task is to solve the wave equation. In axial 
symmetry with II along z everywhere, the wave equation 
becomes: 



1^ 

C2 df^ 



]_d_ 
J.2 Qj. 



n = 0(A13) 



when n, the magnitude of II, has spherical symmetry, 
the solution to the above equation describes a wave emerging 
from the centre and expanding outwards. It is 



where D{ct — r) is any function. 



(A14) 



V X B = E 



(A4) 



The other two equations will let us determine function 
n. Substituting into V • B = we take 

V • (ft - V^n) = 0, (A5) 

and by substituting in the equation V x E = — B we 
have: 



(n- v^n)- =0. 



(A6) 



Given that the fields that have physical meaning are ex- 
pressed in forms of derivatives of II we can choose it to sat- 
isfy the wave equation without overconstraining the fields. 
The proof that follows uses a gauge transformation of II. 



n = n + V X c, 

for C = 0. Indeed: 
E = -V X n = -V X ft = E, 



(A7) 



(A8) 



B = -n + V(V-n) = B. (A9) 
Integrating equations (A5) and (A6), we take: 



APPENDIX B: SELF-SIMILAR MULTIPOLES 

It is possible to extend the idea of the Hertzian dipole to de- 
scribe a multipole. These multipoles have electric and mag- 
netic fields confined between co-axial cones. It is possible to 
isolate and apply the solutions only within the central cone 
or between any two cones ignoring the rest of the field as 
there is no interaction between them. In order to find such 
solutions, the vector potential 11 will be a function of r, t 
and 6. Although it is possible to find general solutions for 
the wave equation in spherical coordinates, we can by de- 
manding self-similarity, design a set of multipole solutions 
that coincide with the Prendergast system. 

By setting II = n(ii, 9)i the wave equation is modified 
to the form: 



d 



(-4) 



n = o, (Bi) 



dv'^ dv sin 9 89 
and by the change of variable it = - it simplifies to 



1) 



d 



sin 9 



n = o. 



(B2) 



When n is independent of 9 the solution is straight- 
forward and equivalent to the case solved in the previous 
section. For ii > 1 we get 
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n = m{u — l)z, 



(B3) 



and the results for the fields are the same as the ones 
described in equation (A13) for a linearly increasing dipolc. 

In the case of a multipole that is for / > it is possible 
to find analytical solutions for equation (B2) by separation 
of variables for m > 1 



Tl = H{u- l)Pi{n){u - l)P/(u)z. 



(B4) 



where H is Heaviside's step function and /u is cos 6. For 
u < 1, zero satisfies equation (B2). In order to find the 
equivalence of these solutions to the Prendorgast fields in 
empty space we will check if the fields defined in (Al) and 
(A2) obey the conditions set in expressions (1) and (6) 



■ V X n = - X 

u 



-n + v(v-n)). 



(B5) 



Indeed (B5) reduces to (B2) after the differentiations 
are done; therefore, a field that obeys the self-similar wave 
equation, it will obey (1) and (6) as well. Thus the Hertzian 
multipole is equivalent to the multipole fields that obey the 
empty space Prendorgast equation. 

To check that this solution (B4) holds not merely for 
M > 1 and M < 1 but also through u = 1 we divide equation 
(B4) by the regular function — (w+ 1) and integrate over the 
small region —e<u — l<e. After integrating by parts the 
left hand side gives: 



/ («- 
Jl-e 



(B6) 



whore H stands for //(m — 1) so it is 1 when u — 1 + e 
and zero when u = 1 — e and its derivatives arc zero at both 
places. The first term reduces to l{l + 1)-Pi(l + e)e and the 
second to — e(2+e)P/(l-|-e) and the third reduces to the same 
integral without the H but with the lower limit replaced by 
1. It is clearly O(e^). Thus the whole expression is 0(e) and 
there is no S function discontinuity in the equation at m = 1. 
It is interesting to remark that had we divided equation (B4) 
by 1 — before integrating we have found the result 



({|^ - i)PUn)] - l{l + 1)[P,(«)]1+')p,(m) (B7) 

which equals + l)Pj(l) = l{l + 1), so the equation 
when divided by 1 — does have a 5 function discontinu- 
ity between the two sides at u = 1! It is only the extra 
factor u — 1 that saves the situation above. The derivation 
of equation (B2) from Maxwell's equations does not involve 
any multiplication hy — 1. Having demonstrated that our 
solution 



n(M, m) = P((m)(w' - l)Pl{u)H(u - 1) (B8) 
is a solution everywhere we recall that our starting 



equations arc independent of the zero point of time so an 
explosion that starts at t — t will be described by setting 
u = c{t — T)/r rather than ct/r. Also, since our equations are 
the linear Maxwell equations, we can superpose solutions to 
obtain, 



/■ro/2 

n(ct,r,M)= / n 

Jro/2 



,c(t — r) .dr 

{- ,w— ■ 

r To 



(B9) 



If H satisfies Maxwell's equations H, will too. The 

physics of this is that the explosion takes a finite time of 
order ro and occurs in a finite volume of order (cro)^ so we 
need to average the fields generated by a point explosion 
over radial distances cto. 

It is not hard to perform this averaging because we are 
interested in regions where u » — with ro small. We 



expand the function (u^ 
u = uo where mo = — • 



l)P/(u) with u 



c(t-T-o) 



about 



(w' - l)P/(w) = (wo - l)f/(Mo) + {u- Uo)l{l + l)Pi{uo) 

+^{u - uofl{l + l)P/(wo) + 0{u - Mo)^(BlO) 

where we used the fact that p satisfies Legendre's equa- 
tion. Now when the step in the Heaviside function in equa- 
tion (B8) is at later times than those involved in the integral 
in (B9) we replace H by unity. So provided r < c(i — to /2) 
we perform the integral in (B9) and find 



H : 



(mo 



„2_2 

c To 



^)PU^o) + + l)PK«o) Pi(/^); (Bll) 



Wo — 1 ^ 



CTo 

2r ■ 



However when \t — r/c\ < to/2 the Heaviside function 

cuts into the range of integration in (B9). If we multiply 
equation (BIO) by H{u — 1) we then find that the range of 
integration has its upper limit changed from to/2 to i — r/c. 
To evaluate fl we need the three integrals 



f-r/, 



ro/2 



(u — Uo) — 
to 



for n=0, 1 and 2. These are readily evaluated: 

' ct — r 



h 



fct-r 1 \ 
^°=(^+2) 

CTo r 1 f ct — r 



2r 

C Tq 

3r2 



CTO 

ct- 

CTo 



(B12) 

(B13) 

(B14) 
(B15) 



Notice that the variable {ct — r)/{cro) lies in the range 
-1/2 to -1-1/2 because \t — r/c\ ^ ro/2 in this region. Thus 



{ul - 1)P/(mo)/o + l{l + l)P(wo)/i+ 



+ ^l{l+'^)Pl{>Mdl2 



Pi{fi),\t-r/c\ ro/2(B16) 



As a further check we took / = 1 and showed that this 
expression for fl does indeed obey equation (A13) which it- 
self follows directly from Maxwell's equations. Thus our pro- 
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cedure of finding singular self-similar solutions and then av- 
eraging out the singularities is indeed giving non-singular so- 
lutions of Maxwell's equations even in the region over which 
the singularity is smoothed out. Of course in the new re- 
gion (III) r ^ c{t + To /2) there is no signal and 11 is zero. 
Remembering that uo = ct/r we see that fl is given by 



^{ul-l) + ^l{l + l)^P;{uo)Pi{n) uo-l>^ 

n = \ {[{ui-i)io + ^Ki + i)i2]p;{uo)+ (B17) 

[ +l{l + l)hPi{uo)}Pi{n) |wo - 1| ^ ^ 

uo-l<-^. 

It is easy to see that fl is continuous as (^7^) = ±1/2 

2 2 

at the edges of region (II), so Jo and I2 go from 1 and 
at the inner edge of region (II), to zero at the outer edge 
and I\ is zero at both edges. Such a continuity is actually 
inevitable since even the unsrnoothod \l(u, ^) is contirmous 
since {v? — l)H(u—l) is continuous but furthermore the first 
derivative of 11 includes a Heaviside function discontinuity. 
The result of smoothing in t ensures that the II is continuous 
but n may have a Heaviside function corresponding to finite 
sudden changes in field strength at the edges of region (II). 
Even such jumps would be avoided had we used the gentler 
function 2^{t-\- ^)'^(t— ^)'^dT with t in the range —to/2 

to -l-To/2 in place of the simpler function ^ used above. The 
more complicated function has zero first derivatives at both 
ends of the range and would lead to continuous and differen- 
tiable field strengths everywhere but wc felt that the much 
more complicated expression for the /„ added unnecessary 
complexity. 

We notice that when ro is small, the smoothing pro- 
duced a negligible modification of order (cTo/r)'^ in region 
(I). In region (II) the changes are significant in that the 5 
function fields are removed and spread out over a region 
Ar = CTo around r = ct. The 7i term is especially signifi- 
cant here, despite its falling to zero at the inner and outer 
edges of this region. 
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